Higgs production via vector-boson fusion at NNLO in QCD 
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We present the total cross sections at next-to-next-to-leading order (NNLO) in the strong coupling 
for Higgs production via weak boson fusion. Our results are obtained via the structure function 
approach, which builds upon the approximate, though very accurate, factorization of the QCD 
corrections between the two quark lines. The theoretical uncertainty on the total cross sections at 
the LHC from higher order corrections and the parton distribution uncertainties are estimated at 
the 2% level each for a wide range of Higgs boson masses. 
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One of the main aims for the now-running LHC col- 
lider machine is to elucidate the mechanism of electro- 
weak symmetry breaking, and in particular to determine 
whether a Standard Model Higgs boson exists or not. 

Among the various production mechanisms for the 
Higgs boson, vector-boson fusion (VBF) offers certainly 
one of the most promising and interesting signals [l|43|. 
The corresponding cross section is second in size out- 
numbered only by the gluon-gluon fusion process, it de- 
creases rather mildly with the Higgs boson mass mn and 
it is proportional to the tree-level Higgs coupling to the 
vector-bosons, W, Z. Moreover, it provides such an ex- 
perimentally clean signature with the presence of at least 
two jets in the forward/backward rapidity region that a 
rich variety of decay modes can be searched for, open- 
ing the access to the very difficult measurements of the 
Higgs couplings Given the importance and the ex- 
perimental prospects for searches of VBF signals it is an 
urgent task to provide the corresponding theory predic- 
tions with the best possible precision including quantum 
corrections. 

Higgs production in VBF is a pure electroweak pro- 
cess at leading order (LO), see Fig. [T] However, at a 
hadron collider the QCD radiative corrections are typi- 
cally sizable, and they have first been computed for the 
total cross-section now almost two decades ago in the so- 
called structure function approach Q . More recently, the 
differential cross-section at next-to-leadingorder (NLO) 
accuracy in QCD has become available [6| along with 
its implementation in a Monte-Carlo event generator Q, 
and also the full set of combined NLO QCD and elec- 
troweak corrections are now known ^. The typical ac- 
curacy of the current QCD predictions can be estimated 
in the 5 — 10% range. 

In this letter we briefly present the results of the 
computation of the dominant contributions to VBF at 
next-to-next-to-leading order (NNLO) accuracy in QCD, 
which give rise to dramatic reductions of the theoretical 
uncertainties. To that end, we are using the structure 



function approach. As we will argue in the following, this 
method, although not truly exact at NNLO, includes the 
bulk of the radiative corrections so that the remaining 
contributions, which are both, parametrically small and 
kineniatically suppressed, can be safely neglected. 

The first point to be addressed is to which extent VBF 
is a well-defined process by itself: interference effects 
with other processes occur, possibly at higher orders in 
the strong and/or electroweak coupling {as resp. asw), 
whose size sets the target accuracy to which VBF as such 
can possibly be determined. For example, VBF processes 
interfere already at LO with the so-called Higgs asso- 
ciated production, leading to a Higgs and two-jet final 
state, i.e., pp — >■ HV'^*^ — >■ Hjj, an effect which is, how- 
ever, at the per mil level Q. At higher orders, interfer- 
ence with gluon-gluon fusion processes can occur starting 
at a\y^a1, but they are also found to be well below the 
percent level 0, [13] • It seems therefore conceivable that 
Higgs production through VBF can be defined in per- 
turbation theory up to an ambiguity of not much better 
than 1%, which sets the target precision for theoretical 
predictions. 

The structure function approach is based on the obser- 
vation that to a very good approximation the VBF pro- 
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FIG. 1: Higgs production via the VBF process. 
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cess can be described as a double deep-inelastic scatter- 
ing process (DIS), see Fig. [U where two (virtual) vector- 
bosons Vi (independently) emitted from the hadronic ini- 
tial states fuse into a Higgs boson. This approximation 
builds on the absence (or smallness) of the QCD inter- 
ference between the two inclusive final states Xi and X2- 



In that case the total cross section is given as a product 
of the matrix element M'^f for VBF, i.e., V/'y/ ^ H, 
and of the DIS hadronic tensor VF^i/, the latter being ex- 
pressed in terms of the standard DIS structure functions 
Fi{x,Q'^) with i = 1,2,3: 
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Here Q| ~ —qf, Xi = Ql/{2Pi ■ qi) are the usual DIS 
variables, Si = [Pi + qiY ^I'e the invariant masses of 
the i-th proton remnant, and denote the vector- 
boson masses, see Fig. [T] Gf is Fermi's constant and 
at a given center-of-mass energy of the collider the 
three-particle phase space is given by the second line in 
Eq. Q. 

The factorization underlying Eq. ([T]) does not hold ex- 
actly already at LO, because interference can occur either 
between identical final state quarks (i.e., uu — ^ Huu) or 
between processes where either a W or a Z can be ex- 
changed (i.e., ud — > Hud). However, at LO, they can 
be easily computed and have been included in our re- 
sults. On the other hand, simple arguments of kinemat- 
ics (based on the behavior of the propagators in the ma- 
trix element [llj) show that such contributions are heav- 
ily suppressed and contribute to the total cross section 
well below the percent level, a fact that has been con- 
firmed by complete calculation even through NLO Q. 
Apart from these interference effects, the factorization of 
Eq. (dJ is still exact at NLO, so that the DIS structure 
functions at NLO can be employed. This is due to 
color conservation: QCD corrections to the upper quark 
line, Figs.[2ja,b), are independent from those of the lower 
line, i.e., Tr(r) = for generators of the color SU(A^c) 
gauge group. 

At NNLO due to the possibility of linking the up- 
per and lower quark lines with two colored particles the 
factorization breaks down. However, it is easy to see 
that such class of corrections. Figs. [3Ua,b), are infrared 
and ultraviolet finite, gauge invariant and suppressed [l3l | 
both kinematically and parametrically by a factor 
The dominant contributions. Figs. [2jc,d) can therefore 
be included in the structure functions at NNLO 



17|. At order a^^ra'^ another class of diagrams arises. 



Figs.|3Uc,d), which contributes significantly to associated 
Higgs production, pp — )■ VH, see While a full 

computation of these diagrams in VBF is in progress and 





FIG. 2: Representative Feynman diagrams of processes in- 
cluded in the structure function approach. 
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FIG. 3: Representative Feynman diagrams of processes not 
included in the structure function approach. 



will be presented elsewhere [20j , a first conservative esti- 
mate can be easily obtained in the rrih — >■ 0,mt -> co 
limit, where the triangle diagrams dominate. In this 
limit, the contribution to the total cross section is well be- 
low the percent level and therefore can be neglected. At 
order a^y^a'^, also di-jet amplitudes with a single quark 
line arise where the Higgs is emitted from a virtual weak 
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FIG. 4: The total cross section at LO, NLO and NNLO as a 
function of mi/ for a ^ = 7 TeV LHC employing the MSTW 
PDF set [23|. The uncertainty bands are obtained by scale 
variation as explained in the text. 
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FIG. 5: The PDF uncertainty of the total cross section at 
NNLO as a function of mn at a ^ = 7 TeV LHC for the 68% 
CL MSTW PDF set For ABKM Q and JR09VF [H 
the ratio of the central value is plotted. 



boson loop. This class of diagrams is gauge invariant but 
not infrared safe and as it is not a VBF process, it is 
not included in our calculation. Its contribution to typi- 
cal VBF final states with a Higgs and two jets has been 



found to be negligible [21 



We now turn to the discussion of the results. For the 
sake of illustration we consider only a V5 = 7 TeV LHC, 
keeping in mind that the conclusions presented here are 
qualitatively the same for a = 14 TeV LHC, and also 
for Tevatron, see (20l |. Our reference parton distribution 
functions (PDFs) set is MSTW Q and the electroweak 
parameters {G p, Mz, Mw, sin'^ 9w) are set to their re- 
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TABLE I: Cross sections (pb) at a ^ = 7 TeV LHC with 
the uncertainty due to independent scale variations fiR,fiF G 
[Q/4,4Q] at LO, NLO and NNLO in QCD as obtained with 
the MSTW PDF sets fH. 



spective PDG values [25 1. 

Fig. m presents the cross section as a function of the 
Higgs mass at LO, NLO and NNLO in QCD, together 
with the uncertainties coming from (uncalculated) higher 
orders. These are estimated by an independent varia- 
tion of the factorization and renormalization scales in 
the range /x_r,/zf = ^r,fQ with ^r^f G [1/4,4], where 
Q is the virtuality of the vector-boson probing the corre- 
sponding structure function to which we apply a technical 
cutoff of 1 GeV. The lower inlay of Fig. |4]zooms in on the 
relative variations normalized to the NNLO cross section 
at ^R,^F = Q, so that the exceptionally good conver- 
gence of the perturbation series can be appreciated. For 
NNLO this is at the 2% level and in principle, could be 
pushed even further within the structure function ap- 
proach by incorporating the available hard corrections at 
order 26l428|. Numbers for our best estimate, i.e., 
NNLO in QCD, are presented in Table 

The most natural choice fJ.R,^J.F = Cr.fQ as a refer- 
ence scale is also supported by kinematics arguments, 
i.e., the observation that the average gauge boson vir- 
tuality in VBF amounts only to (Q) ~ 20 GeV for a 
VS ^ 7 TeV LHC. Of course, other scale choices, e.g. 
IJ'Rt^-f G [?Ti///4,4TO/f], are equally valid. However, they 
typically exhibit a much poorer convergence of the per- 
turbative expansion and lead to sizable deviations in the 
lower order predictions, especially for heavy Higgs bosons 
(e.g. a 7% difference for mjj = 400 GeV at NLO). Only at 
NNLO, both the central values and the uncertainty band 
for the latter choice agree within the 2% level with those 
in Table (U). This clearly demonstrates the markedly im- 
proved scale stability of our NNLO predictions. 

In Fig. [5] the dependence on the parton distribu- 
tions and their errors is studied, which estimates the 
uncertainty of the total cross section due to the non- 
perturbative parton dynamics inside the proton. To this 
aim we employ the MSTW 68% confidence level PDF 
sets 221 through NNLO and compare also with the cen- 
tral predictions obtained with the other available PDF 
sets based on complete NNLO QCD predictions, i.e., 
ABKM [23] and JR09VF [241]. The resuhs are consistent 
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TABLE II: Cross sections (pb) at a = 7 TeV LHC with 
the same parameters as Table ([I]). The approximations NLO^ 
(NNLO^) employ the complete NLO (NNLO) expressions for 
the structure functions both on the upper and lower quark 
lines, cf. Fig. [T] 

and show that an almost constant 2% PDF uncertainty 
at NNLO can be associated to the cross sections for a 
wide range of Higgs boson masses, which is shghtly re- 
duced compared to the NLO value of 2.5%. The PDFs 
are dominantly sampled at an average (x) ~ 0.05 at a 
y/S = l TeV LHC, that is in a region where further non- 
perturbative corrections such as higher-twist effects are 
negligible, see e.g., Ref. 0]. 

Let us next turn to a discussion of the systematics 
of our approach. In Table (jlH) we present an alterna- 
tive study of the perturbative series, namely we consider 
the effects of "improving" the as and expansions of 
the cross section, in terms of expansions of the structure 
functions at NLO and NNLO which enter as a square. 
The NLO^ are obtained by keeping the terms in both 
structure functions, i.e., including formally a1 terms in 
the cross section and using the NLO PDFs, while the 
NNLO^ results perform the same procedure at one higher 
order in ofg . A comparison of the numbers in Tables (U) & 
shows that the different higher order approximations 
agree extremely well within the respective errors bands. 
Implicitly, our findings also demonstrate that the impact 
of the NNLO QCD corrections at the central value, i.e., 
A*ii)A*-F = Q^ is relatively small, and in fact smaller than 
1%, as we have explicitly checked. Altogether this gives 
further evidence that non-factorizable contributions at 
NNLO which are not accounted for in the structure func- 
tion approach and can very conservatively be estimated 
to amount to less than 10% of the NNLO factorizable 
ones, are completely negligible. 

To summarize, we have presented the first computation 
of the VBF cross section at NNLO in QCD. The inclu- 
sion of higher order corrections stabilizes the results at 
the 2% level against arbitrary variations of the renormal- 
ization and factorization scales, indicating an extremely 
well-behaved perturbative expansion. PDF uncertaines 
are estimated at the 2% level as well, uniform over a wide 
range of Higgs boson masses. Our results motivate the 
calculation of the differential cross sections distributions 
at NNLO via an exclusive method such as that proposed 



in Ref. [s^. Moreover, our approach can be used to pro- 
vide cross section predictions to NNLO accuracy for any 
weak-boson fusion process leading to a weakly interact- 
ing n-particle final state X, i.e., V*V'* X. Work in 
the latter direction is in progress. 
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